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Abstract. In nuclear physics, the relativistic mean-field theory describes the 
nucleus as a system of Dirac nucleons which interact via meson fields. In a 
static case and without nonlinear self-coupling of the a meson, the relativistic 
mean-field equations become a system of Dirac equations where the potential 
is given by the meson and photon fields. The aim of this work is to prove the 
existence of solutions of these equations. We consider a minimization problem 
with constraints that involve negative spectral projectors and we apply the 
concentration-compactness lemma to find a minimizer of this problem. We 
show that this minimizer is a solution of the relativistic mean-field equations 
considered. 



1. Introduction 

In this paper, we present the first mathematically rigorous result concerning the 
existence of solutions of the relativistic mean-field equations of the atomic nucleus 
in a static case and without nonlinear self-coupling of the a meson. 

Though often used in practice, the models of nuclear physics have rarely been 
considered from a mathematical point of view: some nonrelativistic models (of 
Hartree-Fock type) were studied by D. Gogny and P.L. Lions in 1986 (PQ), but, to 
our knowledge, there are no rigorous mathematical studies of relativistic models, 
which are however extremely important in nuclear physics. 

In nuclear physics, the relativistic mean- field RMF theory describes the nucleus 
as a system of Dirac nucleons which interact in a relativistic covariant manner via 
meson fields. During the last years, the relativistic mean-field theory has received 
wide attention due to its successful description of lots of nuclear phenomena. The 
relativistic mean-field model is considered to be the relativistic generalization of the 
nonrelativistic models such as the Skyrme force or the Gogny force Hartree-Fock 
theory, using effective mesonic degrees of freedom rather than instantaneous forces. 
The relativistic model describes successfully the single-particle structure of nuclei 
as the nonrelativistic ones and provides a natural explanation of some relativistic 
effects as the spin-orbit force (see [2] , [3] , [4] , [5] ) . 

The model is formulated on the basis of two approximations, the mean-field and 
the no-sea approximation. Thanks to the mean-field approximation, the fields for 
the mesons and the photons are treated as classical fields and the nucleons behave as 
noninteracting particles moving in these mean fields. This implies that the nucleon 
field operator can be expanded in single-particle states ib Q (x^), 



i|> = ^iMx")a Q (1-1) 
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where d a is the annihilation operator for a nucleon in the state a, while the 
densities become simple bilinear sums over the ~\\> a - The no-sea approximation 
corresponds to neglecting the vacuum polarization, that means that we have a 
number of occupied single-particle orbitals il) Q , a = 1, which determines 

the densities. We remind that when the isospir|f| of the particles is not fixed, 
ib Q (x^) G C 2 ® C 4 . Moreover, the single-particle wave functions have to satisfy 
the constraint L 3 lb* (t, x)\^p(t, x) d 3 x = 5 a p- 

The Lagrangian density of the RMF theory can be written as 

^nucleons ~t~ ^mesons ~t~ ^coupling ■ (^■•^J 

The free Lagrangian for the nucleons is 

n 

C nuc leons = ^ W ^a{i (3-2 ® 7^) &a ~ W b )l|> a (1.3) 

a=l 

where m& denotes the nucleon mass, 7 M are the Dirac matrices, w a are occupation 
weights, < w a < 1, and ij) Q = tl>* (l 2 ® 7 ) with 1 2 = '' ' " 



The Lagrangian for the free meson fields is 

1 

2 

1 

2 
1 



1 



--(dWdnUiv-mluVwn) 



1 
2 



—d^A v d l j,A v (1.4) 



where cr, w M and H M describe respectively the a, u> and p meson field, and stands 
for the photon field. Moreover, an antisymmetrized derivative is defined via 



3M" = d^A" - d v A fl . 

We remind that the a meson is an isoscalar scalar meson which provides a medium 
range attractive interaction, the u> meson is an isoscalar vector meson leading to a 
short range repulsive interaction, the p meson is an isovector vector meson needed 
for a better description of isospin-dependent effects in the nuclei, and the photon 
describes the electromagnetic interaction. 
Finally, the Lagrangian for the coupling is 

Coupling = -ga<JPs - 9^ Pa - 9 P R* ■ Pa - eA^ p^ - U (cr) (1.5) 

where U (cr) = ^^c 3 + \b^,a A represents a nonlinear self-coupling of the a meson. 
Note that in Reinhard's paper [2] the coupling constant for the p meson is 2g p . 



a Isospin (contraction of isotopic spin) is a quantum number related to the strong interaction. 
Isospin was introduced by Heisenberg in 1932; he observed that the neutron is almost identical to 
the proton, apart from the fact that it carries no charge. In particular, their masses are close and 
they are indistinguishable under the strong interactions. So, the proton and the neutron appear 
to be two states of the same particle, the nucleon, associated with different isospin projections 

(0.0)- 
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The densities are 



Ps = ^2w a i\) a i\> a , (1.6) 

a=l 

a 

= ^2 w aA>a (12 ®7p)4'a, (1-7) 
a=l 

n 

= X! w q tJj q (t (g> 7^)1^0, (1.8) 

a=l 

= H Wct ^ Q f J( 1 2 + r )(8)7M)^a. (1.9) 



a=l x 7 

We remind that R and p are vectors in isospin space and • denotes the vector 
product therein, and f is the vector of the Pauli matrices which occurs in the 
definition of the isospin operator. More precisely, the three components of the 
isospin operator are defined by t = |r and, in particular, the third component is 
given by 

t 1. 1/10 



2 T ° 2 V I 

I ho pr< ''."D si ate. r< ■])!'< -s< mi 1.1 -! 1 1 >y the vector ( ^ ) . is tlx 1 rij-eiiM.ale of t„ associated 


1 





with the eigenvalue tq — 1 and the neutron state, represented by the vector 



is the eigenstate of f associated with the eigenvalue r = — 1. 

The model contains as free parameters the meson masses m a , m w and m p , as 
well as the coupling constants <? CT , g u , g p , 62 and 63. For the nucleon mass mt> the 
free value is usually employed. 

Most applications of the relativistic mean-field model are concerned with station- 
ary states; then, like in [2], we want to derive the field equations for the static case. 
Moreover, we remark that it is generally true that proton and neutron states do 
not mix, that means that the single-particle states are eingenstates of the operator 
fo. As a consequence, only the components with isospin projection appear, i.e. 
Ron and po M . 

Stationarity implies that all time derivatives and also the spatial components of 
densities and fields vanish; only the fields a, ljq, R 00 and A remain and they are 
independent of time. 

Furthermore, the single-particle wave functions can be written as il) Q = f J 

ipa for protons, and il> Q = ( ^ J ® ip a for neutrons. Each function ip a may be 
separated as 

ll> a (t,x) = e- ie «VaW (1.10) 

where the e a are the single-particle energies and e a > 0. 
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Varying the action integral S = J C d^x with respect to the wave functions and 
to the fields with all the above simplifications inserted yields 



Saloon 

{-A + ml)a + U'{a) 
(-A + m*)wo 
(-A + m 2 p )Roo 
-AAn 



[-i'y ■ V + m b + g a cr + g^ofo 

+9 P RoaloTo + ieA 7o(l + r ) tp a , (1-H) 

-9aPs, (1.12) 

9uPo, (1.13) 

SpPoo, (1.14) 

epg, (1.15) 



with 7 = (7 1 ,7 2 ,7 3 ). This set of equations, together with the definition of the 
densities, constitutes a self-consistent field problem that can be solved numerically 
using an iterative scheme (see [2] , [2] ) • We observe that there is no proof of con- 
vergence of this algorithm. 

In this paper, we consider the case without nonlinear self-coupling of the a meson, 
i.e. &2 = &3 = 0, and we choose a fixed occupation, that means that the occupation 
weights w a are defined as 



■ii'o 



1 a = l,...,A 
otherwise 



(1.16) 



where A is the nucleon number. In this case, the equations (|1. 12111. 15| can be solved 
explicitly and we obtain 



a = — 



9a_ 

4tt 



u 

Roo 
A 

Hence, the equation (|1.11[) becomes 



9lo 
47T 

9p_ 
4tt 

e 

4tt 



*P0 



-rn p \- 



* Poo 



1 



(1.17) 
(1.18) 
(1.19) 
(1.20) 



e a ijj 



+T 



4tt 



A- 



^ s ' + 47 



1 



*Po 



*P00) +^(l+T ) — 



■ Po 



i> a (1.21) 



where Ho = —ia ■ V + firrib is the free Dirac operator, 







f 










) 



for k = 1, 2, 3, with 



1 

1 



02 = 



-i 

1 



^3 



1 

-1 
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The operator Ho acts on 4-spinors, i.e. functions ip E % := L 2 (R 3 ,C 4 ). It is 
self-adjoint on U, with domain if 1 (R 3 , C 4 ) and form-domain E := i? 1 / 2 ^ 3 , C 4 ). 
Moreover, it is defined to ensure 



H 2 = -A + m 2 b . 



The spectrum of Ho is (— oo, —nib] U [m(,,+oo), and the projector associated with 
the negative (resp. positive) part of the spectrum of H will be denoted by A~ 
(resp. A + ). Finally, we endow the space E with the norm \\ip\\ 2 E := (ip, \H \ iP)l 2 - 

Using the convention r = 1 for the protons and t = — 1 for the neutrons, the 
densities can be written as 



Po = 



Poo = 



Po 



fe=i 

A 



fe=i 

z 



k=Z+l 



k=l 
Z 



J2^*kA 



k=i 



(1.22) 
(1.23) 
(1.24) 
(1.25) 



with Z the number of protons, N = A — Z the number of neutrons and ipi = ip*(3; 
furthermore, the nonlinear Dirac equations are given by 



Ho -13 



47T 



9l 

47T 



,-m a \- 



* Pi 



*Poo }+^\—*Po 



9^ 
1 



,-m u \-\ 



*P0 



t/ji = e l ip l (1.26) 



if 1< i < Z, and 



H n ^in := 



q 2 



9 JL 

47T 



* Poo 



9Z / e 



*^' + 47 



Ipi = Silpi 



(1.27) 



if Z + 1 < i < A, with ^ = (?/>i, . . . , ipz,ipz+i, ■ ■ ■ , iPa) and under the constraints 
Jr3 i>ii>j = $ij for 1 < h 3 < Z and for Z + 1 < i,j < A. 

In what follows, V p ^ and V^,* denote the potentials of the nonlinear Dirac equa- 
tions, namely V^,* = H^ — Ho for fi = p,n. 
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Note that the scalars £% can be seen as Lagrange multipliers; indeed, the nonlinear 
Dirac equations are the Euler-Lagrange equations of the energy functional 

= ill WWi-T-f I P f )Ps{ f e-^-y\dxdy 
^J R 3 8?r J J R 3 XR 3 \x-y\ 

+ gi f [ MM e - mu[x - y \ dxdy 

8tt J 7 R 3 XK 3 \x - y\ 



9 2 P f [ Poo(x)Poo(y) 



8tt J 7 R 3 xR 3 \x - y\ 



e -m. P \x-y\ dxdy 



e 



U 



/ / "7 dxdy (1.28) 

8tt J 7 R 3 XR 3 \x - y| 

under the constraints L 3 ift*ipj = &ij for 1 < i,j < Z and for Z + 1 < i,j < A. 
Here we can suppose that the matrix of Lagrange multipliers is diagonal because 
of the fact that £(^) is invariant under the transformations of (ipi, . . . , V>a) of the 
form 

' U p 

U n 

where U p (resp. U n ) is a Z x Z (resp. N x N) unitary matrix. In the energy 
functional, we remark that only the a meson provides an attractive interaction. 
Indeed, if / is a real function, 

MiM e -^-y\ dxdy = C [ l^dfc 
\x - y\ J R 3 k z + V 

with C a positive constant and / the Fourier transform of /. As a consequence, 
the term 

_g± r r eMpM e - m ^- yldxdy 

Stt J Jr3 xR 3 \x-y\ 
is negative and describes an attractive interaction. 

Since the functional (|1.28[) is not bounded from below under the constraints 
J R 3ip*ipj = Sij, as in [5] (see also [S]), we introduce the following minimization 
problem 

I = inf {£(¥);* G (H^ 2 ) A , f tftj = Sy, 1 < i,j < Z, Z + 1 < i,j < A, 

I JR 3 

A-„(^i,...,Vx) = 0, A-„(^ +1i ...,Va)=0} (1.29) 
together with its extension 

I (Ai, . . . , Xa) = inf/^^e^/y, f i>*ip j = \ i 5 ij ,l<i,j<Z, 

I JR 3 

Z + l < ij < A, A-^i,...,^) = 0, 
A-^(Vz+i,...,V'a) = o} (1.30) 

where, for /i = p, n, = X(-oo.o)(^m,*) ^ s the negative spectral projector of the 
operator H^y, 

Ap>(^i,---,^z) = (A~^i,---,A~^z) = A~^* p 
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and 

The idea of using a constraint of the form A~^,\E' Al = 0, for \i = p,n, is due 
to M.J. Esteban and E. Sere in the case of the Dirac-Fock equations (voir [8]). 
This constraint has a physical meaning; more precisely, if we neglect the vacuum 
polarization, the Dirac sea is represented by the negative spectral projector A~^,. 
Indeed, according to Dirac's original ideas, the vacuum is composed of infinitely 
many particles, which completely fill up the negative spectral subspace of H p ^,\ 
these particles form the Dirac sea. So, by Pauli exclusion principle, the single- 
particle energies Si should be strictly positive and, as a consequence, ^ p should 
be in the positive spectral subspace of for [i = p, n. On the one hand, the 

use of the constraint A~^,^> p = is very helpful since it transforms a strongly 
indefinite problem into a minimization problem; on the other hand, dealing with 
this constraint is the main difficulty of the proof of our results. 

In this paper, we prove that, for g a , g^, g p and e sufficiently small, a solution of 
the equations (|1.26j) and (jl.27j) can be obtained as a solution of the minimization 
problem ([Qty . 

Theorem 1.1. If g a , g u , g p and e are sufficiently small, a minimizer of il.29\) is a 
solution of the equations il.26\) and \1.21^ . 

Moreover, the application of the concentration-compactness method ([ID], [TTj ) 
to the minimization problem (|1.29|) yields the following theorem which is our main 
result. 

Theorem 1.2. If g ai g u , g p ande are sufficiently small, any minimizing sequence of 
\1.29\ ) is relatively compact up to a translation if and only if the following condition 
holds 

I < I (Ax, . . . , X A ) + I (1 - Ai, . . . , 1 - X A ) (1.31) 

A 

for all A fe G [0, 1], k = 1, . . . , A, such that J2 G (0, A). 

fc=i 

In particular, if il.31]) holds, there exists a minimum of il.29]) . 

This result is relevant both from mathematical and physical point of view since 
it provides a condition that ensures the existence of a ground state solution of 
the equations (ll.26|) and (jl.27p . Furthermore, this is the first result relating the 
existence of critical points of a strongly indefinite energy functional to strict con- 
centration-compactness inequalities. 

The condition ga,gui,g p and e sufficiently small means that we are in a weakly 
relativistic regime. In our proof of theorems 11.11 and 11.21 this condition is required 
for several reasons. First of all, if g a ,gu),g p and e are sufficiently small, we can 
show that H p <s, is a self-adjoint isomorphism between H 1 / 2 and its dual H^ 1 / 2 , 
whose inverse is bounded independently of VP . Moreover, we need this condition to 
prove that a minimizing sequence of ([L29]) is bounded in {H 1 / 2 ^?)) . We remark 
that the estimates on g&,guj,g p and e are explicit up to this point. Finally, in both 
theorems, we have to apply the implicit function theorem with g a ,{]ui,gp and e as 
parameters. 

This result is different from that obtained by Esteban-Sere on the Dirac-Fock 
equations (see [12], [8]). In [12], by a more sophisticated variational method, 
Esteban-Sere found a infinite sequence of solutions of the Dirac-Fock equations 
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and, in [8], they showed that, in a weakly relativistic regime, the " first " solution 
of the Dirac-Fock equations found in |12j can be viewed as an electronic ground 
state in the sense that it minimizes the Dirac-Fock energy among all electronic 
configurations which are orthogonal to the Dirac sea. Their variational method 
takes advantage of the fact that the Dirac-Fock energy functional is not transla- 
tion invariant: it contains an attractive interaction term, due to the nucleus, which 
confines the electrons. The nonlinear interaction is rather purely repulsive so that 
the use of concentration-compactness is not necessary. On the contrary, the energy 
functional that we consider is invariant under translations and one of the nonlin- 
ear interaction terms is attractive; because of the translation invariance, we are 
naturally led to use the concentration-compactness argument. 

In section [2j we introduce some useful properties of the potential V^,* and of 
the operator for /i = p, n. In section [3l we show how we can apply the 

concentration-compactness argument to the minimization problem (|1.29p . Finally, 
in section [4J we prove theorem 11.11 



2. Properties of the potential V^y 

In this section, we describe some useful properties of the potential and we 
give a condition on the parameters (g a , g u , g pi e, N, Z) which implies that is 
a self-adjoint isomorphism and its inverse is bounded independently of ^ . 

Lemma 2.1. For any 

Vp,* G Z7(R 3 ), 3 < r < oo 
Vn,y G L r (R 3 ), 1 < r < oo. 

Proof. The proof of this lemma is an application of Young's inequality : if / G 
LP(R 3 ), g G L 9 (M 3 ), then 

||/*$||^<||/IWMU« 

with 1 + ± = i + ±, 1 < p, q, r < oo. 

We remark that if * G (i? 1/2 (R 3 )) , then p s , po, poo and pg are in L P (M 3 ) for 
1<P<|. 

Furthermore, using the definition of the Gamma function, we can show that, for 
any A > 0, G L q (R 3 ) for 1 < q < 3. 

Finally, we observe that jij can be written as jij = h\(x) + h,2(x) with hi G L a (M 3 ) 
for 1 < a < 3 and h 2 G L^(M. 3 ) for 3 < (3 < oo, where hi(x) = A for \x\ < 1, 
h\{x) = otherwise. 
Hence, 

2 fp-m a \-\ \ 2 / e -m u \-\ \ Q 2 / -m p \-\ \ 

v »' - -'l(-"H(T*T*(Tr'*) 

+ £ {Kix) *Pt>)+^ (h2(x) * (Q G L r (R 3 ) 

^ V ' ^ v ' 

eLHR 3 ) n irc ( K3 ) GL^R 3 ), 3<r<r c 
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for 3 < r < r c with r c = 9 ~ 3e for any e > 0, and 

- -4 (n^r *"•) + 1 ( f rf*"°) - i (tt**) e r(R3) 

for 1 < r < r c . 

□ 

For reader's convenience, let us remind the following lemma which lists some 
properties of Hq and coulombic potential V(x) = ^j. 

Lemma 2.2 ( 12 J. The coulombic potential V(x) — A satisfies the following 
Hardy-type inequalities: 

(<p,Q**V)<p) ta < \ + |J (v,|Ho|v)z>, (2-1) 

/or aZ/ iy9 G A + (_ff 1 / 2 ) U A _ (iJ 1 ^ 2 ) and for all probability measures fj, onl 3 . More- 
over, 

(<p,(»*V)ip) L z<^(<p,\H \<p) L2 , VcpeH 1 / 2 , (2.2) 

\\(n*V)<p\\ L2 <2\\Vcp\\ L2 , M V EH\ (2.3) 

In the particular case where /i is equal to the Dirac mass at the origin Sq, we 
refer to Burenkov-Evans ([13]) and Tix ([H], [15]) for the inequality (I2.ip . to Herbst 
([16]) and Kato ([17]) for ([221 and to Thaller's book ([18]) for the standard Hardy 
inequality (|2.3|) . The extension of (|2.1|) . (|2.2j) and ()2.3[) to a general probability 
measure /z is immediate. 

Then, using lemma [2~2l and proceeding like in |12j (Lemma 3.1), we obtain the 
following estimates. 



Lemma 2.3. Assume that 

glA + g 2 p max(Z,N) 

4-7T 

glA + glA + g 2 p Z + e 2 Z 

47T 

gjA + gjA + g 2 p N 

47T 



< ^2 + 2^' (2 - 4) 

< _,o 2 o,. > (2-5) 



tt/2 


+ 2/7T 




2 


tt/2 


+ 2/7T 




2 


tt/2 


+ 2/7T 



< "r7o~T~o7Z" (2-6) 



There is a constant h^ > 0, swc/i £/ia£ /or an?/ * G (H^ 2 (R 3 )) such that 

Gram i2 (*) < 1, 



and i> G H 1 ' 2 ^ 3 ) 



hMWm/* < \\H^\\ H ^/2 (2.7) 

with n = p,n. In other words, H^.q, is a self-adjoint isomorphism between H 1 ^ 2 
and its dual H -1 / 2 , whose inverse is bounded independently of "J . 



Finally, a straightforward application of the inequality (|2.3|) yields the following 
lemma. 
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Lemma 2.4. Assume that 



d„ = M±4±«d <1 . (2 . 9) 



For any * G (ff 1/2 (R 3 )) A such that Granijr,j(*) < 1, 



< 4 /2 |i/ | (2.10) 
{1 - d^lHol < \H^\ (2.11) 



for [i = p, n. 



Remark 2.5. Our estimates are far from optimal. In particular, we do not give any 
condition on Too-, to w and m p . We can expect that taking into account the meson 
masses, one can obtain better estimates. 



3. Proof of theorem 11.21 

This theorem is an application of the concentration-compactness argument (see 

[10] . [TTjV Like in (T], if (ip k , . . . is a minimizing sequence of (|1.29[) . then we 

apply the lemma below (proved in |10| ) with the probability Pk in R 3 whose density 
A 

is \p k and p k = ■ 

<=i 

Lemma 3.1. Let {Pk)k be a sequence of probability measures on R . Then there 
exists a subsequence that we still denote by Pk such that one of the following prop- 
erties holds: 

i. (compactness up to a translation) 3y k G l", Ve > 0, 3R < oo 

P k (B(y k ,R)) >l-e; 

ii. (vanishing) Vi? < oo 



sup P k (B(y,R)) ->0; 

y&a N k 

Hi. (dichotomy) 3a G (0,1), Ve > 0, VM < oo, 3R > M, 3y k G R N , 3R k -> 

k 

-hoo such that 

\P k (B(y k ,R )) -a\ < e, \P k (B (y k , R k ) c ) - (1 - a)\ <e. 

In the following subsections, we prove that if the condition (|1.31l) holds, then we 
can rule out dichotomy and vanishing. 

First, we make a few preliminary observations; let ^f k = (ip k , . . . , ip^) be a 
minimizing sequence and g a , g^, g p and e such that d^ < | for /i = p,n, then ty k is 

bounded in (ff^^R 3 )) . Indeed, since ^ is a minimizing sequence, there exists 
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a constant C such that 



A Z 



3 ) L 2 



3=1 



1 A 

9 E (V>*, ^n,**V>. 



3=1 
3 ) L 2 



j=Z+l 



L 2 



3=1 



3=2+1 



1 1 



L 2 



3=1 



3 = 2+1 



Then, using the fact that, for any k S N, t/>* = A + i5fk ipj for 1 < j < Z, ^ 
A+^fcVj for Z + 1 < j < A and the inequalities (12TTU1) and (l2~TTj) . we obtain 



C>£($,|fl- p ,«*h#) £a + ]T (^,|H n ,»*|^ 



3 / L 2 



3=1 



3=2+1 



3=1 3=^+1 

z r ^i/ 2 l A 

3=1 L J 3=^+1 



L 2 



(l-d„) 1/2 



,1/2 



\m\ 2 



As a conclusion, if 2(1 — d^) 1 / 2 — d\l 2 > 0, that means < |, then ||4 ,fe | 
is bounded independently of k and J is bounded from below. 



2 



3.1. Dichotomy does not occur. If dichotomy occurs (case hi.), then, roughly 
speaking, \l/ fc = (ip k , . . . , ip^) can be split into two parts that we denote by = 
(■01 .i, and ^2 = (V 7 * 2i • • • > ^! 2)- More precisely, let £, £ be cut-off func- 

tions: < f '< 1, < C < 1, C(^) = 1 if M < 1, = if \x\ > 2, ((x) = if 



M < 1, C(as) = 1 if M > 2, e,C e P(M 3 ) and let £„, Cm denote £ (-), C (-)• 



We 



set 



8 

= Ci(--y fc )^ fe 

2 

with i?fc — > +00. We remind that dist (supp tp^, supp 7/^2) > 0, 



|#- (0m+^2) I 



L" 



for 2 < p < 3, and 



(see ma, us]). 

Next, we may assume that 



^-(^1+^*2)11^^0 



(3.1) 
(3.2) 

(3.3) 
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A 

for 1 < i,j < Z, Z + 1 < i,j < A, with A* G [0,1] such that J] A l G (0,4). 

i=l 

In fact, suppose that <d k = (9\,...,6 A ) is a minimizing sequence for which the 
dichotomy case occurs. We remind that £ ('J) is invariant under the transformations 
of (^i, . . . , V>a) of the form 



U = 



u p 
o u n 



where U p (resp. U n ) is a Z x Z (resp. x iV) unitary matrix; then, using this 
kind of transformations and writing <d\ = (Spi)0„i), it is clear that we may 
diagonalize Gram^ (9^ J and Gramma (8 pl ). In particular, we have 

Gram La (0* J = diag(A^, . . . , A|), 
Grani^G^) = diag(A| +1 , . . . , X k A ) 

with < A* 5 < 1 for all k G N. Since {A^}fe is a bounded sequence in R, up to a 
subsequence A* 1 — > A, with < A, < 1, and 

k 

Gram^ (e^j.) -Ap.i -> A Pi i = diag(Ai, . . . , A z ), 
Gram i2 (e^ x ) - A n;1 -> A„,i = diag(A z+ i, . . . , Xa)- 

k 

Moreover, as a consequence of the definition of 6 \ and 9 2 , we have 
Gramme* 2 ) - A p , 2 -> A p , 2 = diag(l - A x , . . . , 1 - A z ), 
Gram L 2(6^ 2 ) - A„ j2 A„. 2 = diag(l - Ax+i, ... ,1 - Xa)- 

k 



(3.4) 



(3.5) 



So, to obtain (|3.3I) . we proceed as follows. 

First of all, we define the sets X p = {i G N : 1 < i < Z}, 2 p> i = {i G I p : Aj = 0}, 
2p,2 = {* G 2p : Xi = 1}, 2p,i = I p \ 2 Pj i and X Pi2 = Z p \ 2T p ,2- 

Second, if i G Zp,i, we replace 8^ with ^ x = and, in the same way, if i G Z p ,2, 

we replace 9f 2 with ipf 2 = 0- 
Next, we denote 



©p,i = ( 




r^ik 
PA 


= Gram L 2(6^ 1 ), 


©P,2 = ( 




s~ik 
P,2 


= Gramme* 2 ), 


A P) i = diag 




A P,2 


= diag(l-Ai) ie i p2 



and we remark that A Pi i and A P;2 are invertible matrices. Furthermore, using (I3.4[) 
and (|3.5p . we can write, for k — > +oo, 



G^ = A P>1 + (1), G* 2 = A p , 2 +o(l); 
then Gp 1 and G p 2 are invertible matrices and 

(G^r 1/2 = A" 1 / 2 + (1), (G* a )-V2 = A" 2 /2 + o(l). 
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To conclude, it is enough to consider a small perturbation of 8^ and 6^2- More 
precisely, we take 



By a straightforward calculation, we obtain 

Gram i 2(*^ 1 ) = A P)1 and Gram L 2 2 ) = A p . 2 ; 
hence, writing ^> k p l = (tp^, . . ., and = (^ 2 , . . . ,V>| i2 )> we have 

Gram L2 (*p x ) = diag(Ai, ...,X Z ) and Gram L 2(*p 2 ) = diag(l - Ai, ... ,1 - X z ). 

Finally, with the same arguments, we can construct = C^z+i n * • • ' l) 
and 2 = (^1+1,2) • - • > ^a,2) such that Gram i 2 (*n,i) = diag(A z+ i, . . . , X A ) and 
Gram^^* 2 ) = diag(l - X z +i, . . . , 1 - X A )- 

We remark that *f = (^fj, . . . , -0a, 1) an< ^ *2 = W*i 2> • ■ • > 2) do not necessar- 
ily satisfy the constraints of / (Ai, . . . , Aa) and I (1 — Ai, . . . , 1 — Aa) respectively, 
then we proceed as follows. 
First, we show that 

A- (^ )1 ,.--^i, 1 W0in (^V2 (M 3 } ^ 
A ,;^(^ + i,i"--^l 1 )->0in (^V2 (R 3 } y 



11 in [ // 1 -it:" T 

fe 

and 



(3.6) 



A-^(^ 2 ,...^i,2)^0in (H 1 ^ 
a;^ (V£+i, 2 , • • • , < 2 ) T> in (HWQP)) 1 



(3.7) 



> U 111 {H*'*{K°)' 

k 

Second, using the implicit function theorem, we construct 

3? = (^,i,<i),^ = (*£ 2 ,< 2 ) e (h 1/2 (m 3 )) Z x (if 1 /2 (R s ) j iV ! 

small perturbations of * 2 in (# 1/2 (R 3 ))' 4 , such that 

A ;,*^p,2 = 0, A- 45 $* i2 = (3.9) 

and 

Gram i2 (3* 4 ) = Gram^ (3.10) 

for /i — p,n and i = 1,2. 

In conclusion, thanks to the continuity of £ in JJ 1//2 (K 3 ), we obtain 
I = lim £(# fc ) > Km £ + Km 

k— >oo fc— >oo 

= lim £($£) + lim 

fc— >oo k^-oo 

> 7(Ai,...,Aa)+/(1-Ai,...,1-A a ) 

that clearly contradicts ()1.31[) . We remind that the first inequality is obtained by 
using the properties of localization of *J , 4^, V^f and Vvt^. 



14 



SIMONA ROTA NODARI 



We start by showing that 

A"^ « l5 . . . , V>|,i) ^ in (HV*(R 3 j) Z . 
Using the formula (see |17j ) 

1 f + °° 

A B - A- = —J [(A - in)- 1 - (B — in)- 1 ] dr, 
1 f +oc 

= 2^J - irf) (B - A)(B - in) dr], 



(3.11) 



we can write 



A P) *^i-A piWfc ^i 



1 

2^ 



{H p ^k - irf) l {H p ^k - Hp ^^iHp ^k - ir)) Vi.i dr}. (3.12) 



for i = 1, . . . , Z. Hence, if we prove that 



and 



i„/.fc 
i,i 



dn — > 0, 



-> 0. 

ffl/2 fc 



HV2 



(Hp^k—ir}) 1 {H p ^ s ,k - Hp^^iHp^k - irf) 1 ip, 
we can conclude that 

First of all, we consider 

f k (v) = (-Hp,** - i^iHp^ - H p ^ k ){Hp^ k - ir})- 1 ^ 
and we prove that, \/rj £ R, 

/ fe (r?) -> 0. 

We decompose the proof of this fact into two lemmas. 

Lemma 3.2. Assume that g a ,gLj,gp and e are sufficiently small and let ^> k be a 
sequence in (H 1/2 (R 3 )) A such that Gram £ 2(* fc ) < 1. Vfc G N, and \\^ k \\ {LP j A is 
bounded independently of k for 2 < p < 3. Then for any (p S L 2 (R 3 ) and for any 
r\ € K, there exists a constant h p such that 



I {Hp ^k — ir/) VI 



< 



1 



IMI 



c 



(h* + ry 2 ) 1/3 



lH1/2 ~ (ml+rf) 1 / 4 
with C a constant that does not depend on k. 
Proof. First of all, we write 

{H p>9 k - ir?)"V = x 

V = (Hp^k - irj)x 

{H a - in)- 1 (Lp - Vp^x) = X 



IMI 



L- 



(3.13) 
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where V p ^k = H p yk—Ho. It is easy to show that if tp 6 L 2 (R 3 ), thenx 6 H 1 / 2 (R 3 ); 
indeed, there exists a constant h p > such that 

\\<p\\l» = (( H p,v -iv)x, (Hp,** -iv)x) L * = \\H p ^kx\\ L 2 + V 2 \\x\\% 
> m b h 2 p \\ X f m/2 + r? \\x\\l, > mlhl \\x\\l, + V 2 \\x\\» 



thanks to Sobolev embeddings and lemma [ 
Next, to have a good estimate of the if 1//2 -norm, we use its definition and we obtain 



\\xf H i/2 



\(H -ir)) 1 (v>-V p ^kX)\\ H i / 2 



< 



(ml + \ P \ 2 ) 1/2 

3 

K + H 2 ) 1/2 

» Ho(p) 2 +V 2 
{ml + \p?) 1/2 



0(p) - Vp^kxip) 



Ho(p) - ir\ 



dp 



<i>{p)\ + v Pi9 hx(p) 



dp 



dp 



< 



< 



m 2 + | P | 2 +?? 2 V 1 ^ 1 + Vp ^ X{P) 

K + H 2 ) 1 / 2 / 2 - — - 
K + b|2) i /2K + 7?2) i /2 \\<P<P)\ + v p ^x(p) 



dp 



(m 2 b +n 2 ) 1 / 2 



2 i |W/ II 2 

L 2 + H^W^X 



To conclude, we have to find an estimate for ||^, i *fcx|| i 2- I n particular, we have 

||Vp,* fc x|| L 2 < H^p.^lliis llxllx,»/4 < ll^p,* fc ILi s IMIi/ [belli/ 



< c I 



1/3 \\ m ., _^-l,J| 2 /3 < 



L 2 



(F**-ir?) d 



C 



L 2 



(/ ) , 2 +7/ 2 ) 1 /3 



IMI 



L 2 



where ft. p = m},h p and C is a constant that does not depend on k. Hence, 

1 f„ , C 



< 



(ml + n 2 ) 1 ^ 



IMI 



L 2 



(/ l 2 + 7 7 2 )V3 



L 2 



□ 



Lemma 3.3. Assume that g a ,gui,gp and e are sufficiently small; let fy k be a mil 
imizing sequence of \1.29\) and Vf^ defined as above. Then, for i = 1, . . . , Z , 



(H Pi n,k - H Pi yk)(H p> v* - iff) Vry G 

Proof. First of all, we study the behavior of 

in R 3 \ B (2/ fc , Writing ^ *(•) = £**(•- we obtain 



(3.14) 



(3.15) 
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with r fc defined by 



(3.16) 



As a consequence, we have 



L- 



( B 3^B(>,%)) 



< 



-fell 



£2(h3^ B (»*,%)) 



< C 



llL 2 (R3) 



Il 2 (r 3 ) 



(3.17) 



since supp ^ yfc = S (y fe , -^). Then, to prove that (H p ^k — irj) converges 
to in L 2 (R 3 \ B (y k , ^j 1 )), it is enough to show that the norm of the commutator 



(H mk - irj) 1 ,^R i 



converges to 0. We remark that 



Hence, using lemma [3~2l we obtain 



l T HffV2(R3) = 0(i? fc 1 ). 



(3.18) 



Finally, using the fact that — > +00, 



(Bp,** - *»?) W,i 



(3.19) 



and 



L*(r=Kb(>,^)) it 



(3.20) 



for 2 < p < 3 thanks to interpolation inequality and Sobolev embeddings. Indeed, 
we remind that (H pi pk — £77)"" is bounded in H 1 / 2 { 
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Second, we consider the potential 

W := Hp-yk — Hp^/k — 



A ' p -m a \-\ 



3 = 1 



a 2 A / e - m "H 
i=i 

.2 ^ 



(w.i a -wr) 



47T ^ 



2 A 



J = Z+1 



e 

47T 



£(o*(W.l'-Wf)) 

.7 = 1 V ' 7 



(3.21) 



and we estimate the L -norm of W in Bf. := B(y , ^-). Using (|3.1|) and the 
definitions of V^i and V'? 2> we obtain, for 1 < p < I, 



J,i " " 
5*1/)* II 



and, in the same way, 



.Ml 



c l 2 



< 



LP(B k 



+ 6-11^ + ^2-^" 



LP(B k 



') k 



Next, we remark that this potential contains three types of terms; for the first one, 
we have 



L 7 (B k ) 



< 



L35/12 (Bfc) 

Similarly, for the second type of terms, we obtain 



m,^i-m)\\ L ^B k) t°- 



L 7 (B k 



< 



-m u |-| 



(w..i 2 -wf) 



L35/12( Bfc ) 



L5/4( Bfe ) fc 



0. 



For the last term, we remind that jij can be written as — hi{x) + fi2(x) with 
hi e L 35/12 (]R 3 ) and h 2 E L 7 (R 3 ), where hi(x) = A for \x\ < 1, hi(x) = 
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otherwise. Then 



c l 2 



L 7 (B k 



< \\h 



l|| L 35/12 (Bfc) 



-> 0. 

Li(B k ) k 



Finally, 



(3.22) 



In conclusion, using (|3.20p and (13.221) . we obtain 

2 

L 2 (R 3 ) 

2 

< 



(Hp^fc - H p -q,k)(H p ^ k - irj) l tp it 



2 

L 2 (R 3 ) 



W*(^ Pl *}-iT?)-Xl 



L 2 ( B 



< \\w k 



W K 



\L'{ B 



Z, 2 (r3 x b(V 

(■Hp,*} - 



Il 7 (r 3 ) 



L"/5(R3) 

2 

L 14 / 5 (R 3 \B(iy fc ,-^)) 



□ 



Hence, if we apply lemma l3~2l to ip = {H p ^h — H p ^k)(H p<5/ k — irj) V%i an d we 
use the result of lemma 13.31 we can conclude that 



{Hp,-** - irj) 1 (H p ^k- H p ^k){H p ^u - irj) V*i 



-> 

H 1 / 2 k 



(3.23) 



for all 77 G M. 

Finally, to prove that f k {rj) d-q — > as k — > oo, we use the Lebesgue's 
dominated convergence theorem. Indeed, the sequence f k converges to / = for 
all 7] £ R and is dominated by an integrable function g. In particular, using lemma 
and its proof, we remark that, Vfc £ N, 



\f k (v)\ < 



1 



c 



< 



(mf+ rfy/ty (^2+^2)1/3 

L + 

(mg + ry 2 )V4(/i2 + ^2)1/3 ^ + ^y/a ^ 



W k (H p ^-ir,)- 1 rp k 1 
C 



L- 



< 



c 



c 



{ml+rfy/^hl + rj 2 ) 1 /* \ {h? p + rf) l l^ j 
and g(rj) £ L X (R). Then 



f k (r))dr,->0. 

k 



W,i\\ L 2 -sW 
= 5(»7) 
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Now, to prove that = A^^Rj. (• - y k )ip k converges to in i? 1 / 2 ^ 3 ), 

we give an estimate on the commutator A~ ^,£.5^ (• — y k ) . Writing k (-) = 

L 8 J 8 

$R k (■ — y k ) and using Cauchy's formula, we infer 

8 

,+oc x 1 



2TT J_, 

%_ f +oc 1 

27T J_ oc Hp -^k 



i r ■ £s±. „k £ R k ..k JT . ■ drj 

H p q,k + it] s >v s >y Hp ^k + irj 



oc #p,*<= +111 V s -V 
+oo t 

——a • V£rj. fc — — drj. 

+ irj s <y Hp tyk + it] 



1 



H p q,k + irj 



di] 



Hence, 



< C 



(to 2 + v 2 ) 1/i {m 2 b + r] 2 ) 1/2 



dri = O(R^). 



Then, as Rk — > +oo, we obtain 
fe 



H 1 / 2 fe 



and, since A p ^i/^i 
conclude that 



A 



p,**)^,^ ^ + £«t ^A^^f and A p>9k ip!? = 0, we 



A " »(^,i,...,^,i)t>0 
p<^i fe 

in (i/ 1/,2 (]R 3 )) Z . Moreover, with the same arguments used above, we prove that 



in (H 1 / 2 (W 3 )) N . 

Furthermore, to show that 



A" 



vfc (^i,2.---.^z,2)T>0 and A niJ/fc (i/.| +12 ,...,Vi 2 ) -> 



p,WJ 

, z 



in (_ff 1 / 2 (M 3 )) z ' and (i/ 1 / 2 (R 3 )) ^ respectively, we can proceed as before; only the 
proof of 

(3.24) 



L 



is slightly different. In this case, 

thanks to the localization property of V^i> and (H p ^k —irj)^ 1 ^^ converges strongly 
to zero in L p (£? (y fe , -^ k )) for 2 < p < 3. In conclusion, 

2 

(■ffp,** - H p,* k )( H p,y* - «?)~ ^,2 L2(R3) 

2 

(-^p,*£ ~~ Hp-^k) 



+ c 2 



l 7 (rkb(/ i)j "fe^ °' 
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Now, we want to construct $f = ($^1, $£,1), *2 = (*p,2' $£,2) e (7f 1/2 ( 
(T/ 1 / 2 ^ 3 ))^, small perturbations of in (Tf 1 / 2 ^))^ 4 , that sat i s f y the con- 

straints of 7 (Ai, . . . , Xa) and 7 (1 — Ai, . . . , 1 — A^) respectively. For this purpose, 
we use the following lemma and its corollary. The proofs of the lemma and the 
corollary are given in the appendix. 



N 



Lemma 3.4. Take * = (* p ,*„) € (Tf 1 / 2 ^ 3 )) x (Tf 1 / 2 ^ 3 ))" sucA tfmi 

z. Gram^ (\J/ p ) := G p < lz and Gram L 2 :— G n < In are invertible 

matrices; 



(jyi/2)Z 



< 5 



( ff l/2)JV 



< (5 



/or <5 > sma/Z enough. 



If g a , gui,g P and e are sufficiently small, there exists $ = (<f> p , $„) e (77 1 / 2 (R 3 )) Z x 



(77 1 /2(R3^ JV such that 



Moreover, 



K.^n = 0. 



Gram L 2 ($ p ) = G p , 
Gram L 2 ($„) = G„. 



(3.25) 
(3.26) 

(3.27) 
(3.28) 



A' 



a sequence of 



Corollary 3.5. Take = (**,**) G (/f 1 / 2 ^ 3 )) x (i7 1 / 2 ( 

functions bounded in (i7 1 / 2 (M 3 )) A smc/i i/i<rf 

z. Gram L 2 (* p ) := G p < \z and Gram £ 2 := G n < In are invertible 

matrices that do not depend on k for any k G N; 



P ,* fc p I 



(ffl/2)Z fe 



0. 



|(JJ1/2)JV fe 

7/ ffcr , 5w , .9,9 a^rf e are sufficiently small, there is a constant fc € N sucft i/iai, /or 
any fe > fc , i/iere exists <& fc = G (i7 1 / 2 ( 

following properties: 

a. 



Z x (7^/2 ( M 3))^ with the 



A * fc $„ = 
p.'f' 1 P 

A" 

n,* fc " 



0, 

= 0. 



b. 



^pii(ffi/2)z y ' 



I <J> fe — V|/ fe II — 5-0 



(3.29) 
(3.30) 

(3.31) 
(3.32) 
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Gram L 2 ($*) = G p , (3.33) 
Gram L 2 = G„. (3.34) 

So, using the corollary |3.5[ we can conclude that if g a ,g^,g p and e are sufficiently 
small, there is a constant fco € N such that, for any k > ko, there exists ^ = 
e (ff 1/2 (R 3 )) Z x (H 1 / 2 ^ 3 ))* with the following properties: 

i. 



A" 



A — rF.fc 



71, $ 



$ te , = 

^n,l w 



(3.35) 
(3.36) 



I cb k _ II — \ n 

| p,l p,l II (i/ 1 / 2 ) 2 77 ' 

I <j) fe _ \i/' c || — > n 

I n,l n,l || (ffi/2\N 77 



iii. For 1 < i, j < Z , Z + 1 < i, j < A, 



XiSij . 



(3.37) 
(3.38) 

(3.39) 



In particular, if 

Gram L 2 (*p 4 ) = diag(A x , . . . , X z ) and Gram L 2 (^^x) = diag(Az+i, . . . , X A ) 
are invertible matrices, we apply the corollary |3.5l to 

On the other hand, if Gram L 2 -J or Gram L 2 (\&„ x j is not an invertible matrix; 
then there exists i <E {1, . . . , A} such that \ — 0. As a consequence, ip^ = for 
any fc £ N. 

We assume, without loss of generality, that A^ = for 1 < i < r p , Z + 1 < i < r n 
and Xi 7^ for r p < i < Z, r n < i < A, and we denote = (tp* 1; . . . , i\) k z x ) and 

Ml = «,i,..., V>Xi). Since 



Gram L 2 = diag(A rp , . . . , A z ) and Gram L 2 = diag(A r 



,Aa) 



are invertible matrices, we can apply the corollary 13.51 to = J J to 

obtain, for any fc > fc , = (&p,i> $n,i) sucn tnat 
i. 



->0, 

( ff l/2)Z-r p + l fc 



n.l 



->0, 

( ff l/2)N-r„ + l fc 



-n,l 

iii. Gram L 2 ($ Pt i) = Gram L 2 and Gram £2 = Gram L 2 
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To conclude, it is enough to take 

**i = (o,...,o,^ il ,...,^| il ) 

and 

and remark that H^ ^k = $ k for fi — p, n. 

In the same way, if ga,guj,gp and e are sufficiently small, there is a constant 
k E N such that, for any k > k , there exists = ($p ;2 , $n, 2 ) e (-ff 1/2 (R 3 )) Z x 
(i/ 1 / 2 (]R 3 )) N with the following properties: 



1. 



A ;* 2 * $ p,2 = 0, (3.40) 



Ar.»< 2 =0. (3.41) 



|*p,a-*p,2|| (1 ,v2)*^ . (3-42) 



-n,2 *ra,2||( ff i/2)N fc 

hi. For 1 < i, j < Z, Z + 1 < i, j < A, 

Ca^i.2 = (1 - Ai)%. (3.44) 

Using (1337)) . (|3~38| . (|3~42| . (|3T43|) and the continuity of £ , we remark that 
lim = lim £(®i), 

k— ¥oo »oo 

lim = lim 



and then, if dichotomy occurs, we have 



/ = lim £(^ fe ) > lim £(^\) + lim £(V h 2 



— , _ ^) 

> °° fc— >0O fc— >0O 



= lim + lim 

fc— J-OO k^-oo 

> I(\ 1 ,...,\a)+I(1-\i,...,1-\a). (3.45) 
It is now clear that (|3. 451) contradicts (II .3 1[) . 

3.2. Vanishing does not occur. If vanishing occurs (case ii.), then Vi? < oo 



sup 

yem.3 JB(y,R) 



for j = 1, . . . , A and ipi, . . . , converge strongly in L P (R 3 ) to for 2 < p < 3 (see 
lemma 7.2 of [ID]). As a consequence, 



3=1 

and 



lim £(* fe ) = V lim / $*H ifi, 

3=1 

A 

/(Ai, . . . , Xa) = m b y~] Xj 



3 = 1 
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thanks to the constraints of the problem. 
This contradicts (|1 .31|) because we have 



A 



I = m b A = m b Aj + mj } y (l - \j) = I (Ai, 

3=1 



A 

.£< 

j=i 



,Aa)+/(1-Ai 



1-A a ) 



At this point, we have shown that any minimizing sequence satisfies the following 
compactness criterion: 3y k 6 R 3 , Ve > 0, 3R < oo 



1 A 

3=1 



> 1 - e. 



B(y k ,R) 



We denote ^ = ^ r * : (- + and we remark that the energy functional £ is invariant 
by translations and $ is in the minimizing set; then $ is a minimizing sequence of 
(TQ9| , Since ^ fc is bounded in (tf 1 / 2 ^ 3 )) A , § ft converges weakly in (H^ 2 (R 3 )) A , 
almost everywhere on R 3 and in (Lj' oc (lR 3 ))' 4 for 2 < p < 3 to some ^; moreover, 
thanks to the concentration-compactness argument, iff k converges strongly to ^> in 
(L 2 (M 3 )) A and in (LP(R 3 )) A for 2 < p < 3. 
As 



for fc — > +oo, it is clear that 



lim 

k— » + oo 



for 1 < i,j < Z and Z + 1 < i,j < A. Furthermore, A ~ W M = for [i = p,n. 
Indeed, as before, 

A~ - A~ - .ibj 



1 

2^ 



and 







since — — > for 2 < p < 3. Then, applying the Lebesgue's dominated 

J k 

convergence theorem as above, we obtain 

A" T fa 

for n=pi{l<j<Z and fj, = ni£Z+l<j<A. As a consequence, 



-> 

ifl/2 fc 



A „ ,f>3 - A „ &H*I>3 



-> 

L 2 k 



and 



A 



L- 



lim 

k— >+oo 



L 2 



lim 

k— >-t-oo 



L 2 



thanks to the properties of the spectral projection A ^ fc and using the fact that 
\\ipi — tfj k \\l2 — >■ 0. So we can conclude that \E" satisfies the constraints of the 

J k 

minimization problem (| 1 . 29[) . 
Finally, we have to prove that 

£(#) < liminf £(* fc ). 

fc— J' + OG 
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It is clear that if Mj - ?/C|| LP — > for 2 < p < 3, then 



(3.46) 



for fj, = p if 1 < j < Z and fj, = n if Z + 1 < j < A. Moreover, we observe that 



< 



(A--A--J(^-^) 



A :„)v^ 



and, with the same arguments used above, we obtain 



A ^ 



= lim 

H 1 / 2 fe->+tx 



A"^ 



(3.47) 



Then, using (|3.46|) . (|3.47|) and the weak lower semicontinuity of the iJ 1 / 2 -norm, we 
get 

A A 

£(*) = ^(^,|-ffo|^) i2 -2^(A-^,|ifo|A-^ 



3=1 



3 = 1 



L 2 



3=1 

< liminf £(^ k ) < £(#). 

k—>-\-oo 



3=Z+1 



L 2 



As a conclusion, W is a minimizer of (|1.29[) and the minimizing sequence ^ is 
relatively compact in (H 1 / 2 ) A up to a translation. 

3.3. The subadditivity condition. To conclude the proof of theorem II. 2 [ it re- 
mains to show that the strict subadditivity condition (|1.31l) is a necessary condition 
for the compactness of all minimizing sequences (see [10], [11]). 
First of all, we prove that we always have 

J</(A 1 ,...,A A ) + /(1-A 1 ,...,1-A A ) (3.48) 

A 

for all A fe € [0, 1], k = 1, . . . , A, such that J2 x k € (0, A). 

k=l 

Let e > and *f , *| be satisfy 

f / (Ai, . . . , X A ) < S (tff ) < I (Ai, . . . , X A ) + e, 
Gram L2 (^ 1 ) = diag (A x , . . . , X z ) , 
Gram L2 (*^ 1 ) = diag (A z+ i, . . . , X A ) , 

and 

/(1-Ai,...,l-A A ) <£ </(1-Ai,...,1-Aa)+e, 
Gram L2 (^ 2 ) = diag (1 - A X) . . . , 1 - X z ) , 
Grainy (1<^ 2 ) = diag(l - X z+ i, ...,1-X A ), 
l A ^|*P,2 = 0, A^.¥» i2 = 0. 

By a density argument, we may assume that "ff and ^ have compact support and 
we denote by = + ^v) where t] is some given unit vector in R 3 . Since for 



(3.49) 



(3.50) 



THE RELATIVISTIC MEAN-FIELD EQUATIONS OF THE ATOMIC NUCLEUS 



2 r , 



k large enough the distance between the supports of and \?|' is strictly positive 
and goes to +00 as k goes to +00, we deduce 



£(^> k ) - £ (*f ) + £ (*l' K ) 



E,k\ 



ipi' ipj -toy 



A 



e.k 



e.k 



— 5" 

-> 

k 



1 < i,j < Z, Z + 1 < i,j < A 
j = 1,...,Z, 



j = Z+l,...,A 



with = *l+%' k - Indeed, A"^,^'" = A-^^-A-^^+A"^^ 



e,fc 



(3.51) 



/ £ >k 



A ye.kipj'^ & n <i, by arguments similar to those used above, we obtain 



k 



0. 



if 1 / 2 fc 







for fj,=p,n. Then, as before, we can construct $ £ ' fc , small perturbation of ^ e,k in 
(if 1 / 2 ^ 3 ))" 4 , such that 



f /(Ai,...,Aa)+/(1-Ai,...,1-Aa) < I™ £ ($ 

fc— >+oo 



e,fc\ 



lim £f* e ' fc ) 



= £ (*f ) + £ (%) < I (A x , . . . , Aa) + /(I - Ax, . . . , 1 - X A ) + 2e, 



Gram i2 ($^) = l z , Gram L2 ($^ fc ) = 1 



(3.52) 



and, by definition of /, we conclude 

I < I (Ai, . . . , Aa) + / (1 - Ai, . . . , 1 - Aa) + 2e. 
In fact, this argument prove also that if 

1 = J(A 1 ,...,Aa) + /(1- Ai,...,1-Aa), 

then there exists a minimizing sequence that is not relatively compact. Indeed, 
let it* and ^2 b e minimizing sequences of I (Ai, . . . , Aa) and I (1 — Ai, . . . , 1 — Aa) 

respectively, with compact support and such that dist (supp tpf L , supp ^2) — * +°° 

' ' fc 

for i = 1, . . . , A. If we take i/jf — V^i+W^) we can show that J R3 ip!?(x)x(x) da; — ^ 

for all x S 2?(R 3 ); then ^ converges weakly to in (_ff 1 / 2 ) yl . Now, as before, we 
can construct $ fe , small perturbation of ^ k in (ff 1 / 2 ) j4 , which is in the minimizing 
set of / and such that 

lim £ ($ k ) =/(A 1 ,...,Aa)+/(1-A 1 ,...,1-Aa) = L 

k^+00 

As a conclusion, is a minimizing sequence that cannot be relatively compact. 
4. Solutions of the relativistic mean-field equations 



In this section, we prove that, in a weakly relativistic regime, a minimizer of 
(|1.29|) is a solution of the equations (|1.26j) and (|1.27j) . 
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X 



Let 

{7 G B(K); 7 = 7*, (ml - A) 1 ^^ - A) 1 / 4 e oi(W)} (4.1) 

where B(H) is the space of bounded linear maps from H to H and a 1(H) is the 
space of trace-class operators on H. 
Now, to each P £ N, we associate 

rp-{ 7 eX; 7 2 =7,tr( 7 ) = P}. 

Given 7 = (7^, 7„) e I x I, wc define 



h °-4 



47T 



Ho-/? 



-m p |- 



9%_ 

47T 



,—m a \ 









47T 




















47T 



-m u |-| 



*P0 



(4.2) 
(4.3) 



7p 

mj-| 



*P0 



(4.4) 



47T 



* Poo 



where 



p s (a;) = p p (x) + p n (x) 

po(x) = p p (x)+p n (x) 

Poo(x) = p p (x)-p n (x) 

with p p {x) = tr(/3 7p (x,x)), p n {x) = tr(^7„(x, x)), p p (x) = ti(^ p (x,x)) and p n (x) = 
tr(7„(a;, x)). 
Finally, we define 

Ap, 7 = Xm±(H P ,-y)> 

Ki,-y = XK±(^n, 7 )- 

Let "J = (^p^n) be a minimizer of the problem (|1.29l) ; to prove that ipi is a 
solution of (H~2%|) for 1 < i < Z and of (TOT)) for Z + 1 < i < A, we proceed as 
follow: first, we consider j p and 7„ the orthogonal projectors defined by 

z 

7p = £|&M&I (4-5) 



and 



.1 



(4.6) 



i=Z+l 

and we denote 7 = (7 P ,7 n ); then, we show that 

[H^,%] = 

for fi = p,n. This implies 

Hp ^i = Si^i for 1 < i < Z, 
H n = £ i^i for Z + 1 < £ < A 
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First of all, we observe that if 'J is a minimizer of (|1.29[) . then the vector 7 
(7p)7n) is a minimizer of the energy 



tr(fr 7p)+tr(fo7n)-|j 



Ps(x)p s (y) m<J \x-y\ 



9u_ 
"8tt 

'8tt 



on N = Y\ x where 



Po( x )Po(v) - 
\x - y\ 

Poo(x)poo(y) < 
N - 2/1 

Pp(x)p P (y) 



; r3 |a: — 2/| 
"' '• y dxdy 



dxdy 



-m p \x—y\ 



dxdy 



<R3 \x-y\ 



dxdy 



(4.7) 



r z = {lp G T z ; 7 P = A+ 7 7 P A+ 7 } , 
r« = {in e r N ;<y n = A+ 7 7„A+ T } 



Next, we remind that 



with H+ n = A+ 7 if jUi7 A+ 7 and = A^H^A^ for fi=p,n. Then 



It is clear that T% N is a subset of 



Fz,n = {7 = (7 P! 7n) G T z x Tat; [i? P)7 ,7p] = 0, [H n ^,i n ] = 0} 



and, since 7 <E N , we obtain £f -, 7^ 



for = p, n. Thus, to conclude, we 



have to prove that H^-, 7^ =0 for p = p,n. We proceed by contradiction. 



We suppose that 



and 



^n~,7 ' 7" 



are different from zero and we define 



7p = K% (K) 1 : = ex P 

7« = Kin (KV 1 : = ex P [ H n^>ln 

In particular, 



#jt 7 . % )% ex P ( e #^7 . % 



In exp e 



»=i 

A 



(4.8) 
(4.9) 

(4.10) 
(4.11) 



i=Z+l 



with ^| = Upipi for 1 < i < Z and = U^ipi for Z + 1 < i < A. 

Using lemma EOI we construct j £ = (7^,7^), small perturbation of i £ = (lp,l n ) 

such that 

<y e = A + 'y 6 A+ 

lp p,7 e "p p,7 e ' 

7 e = A+ e 7 e A + g 
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We remark that 7 E 6 and 



!>?}< 

1=1 

E 10? 



where 



with $f 

A- 

that 



with B u 



i=Z+l 

l...,^ z ) = % = % + + %- + 0(e 2 ) 
. . . , - = + + $r + 0(e 2 ) 

1/2 



(4.12) 
(4.13) 



(4.14) 
(4.15) 



Gram L2 (A+- E ^) 



for /i = p,n. Finally, we remind 

(4.16) 



1/2 



Grainy (A + - e ^f e ) for fi = p,n (see the proof of lemma [3T4] 



Then, to show that we have a contradiction, we want to prove that 

£{l £ p,ln) < £(7j»7n)- 

For this purpose, we calculate £(7^,7^) — £(7 P ,7„); since (7^,7^) is a small per- 
turbation of (7 P ,7n), we can write 

£(%,ln) ~ £(%, 7») = tr (i? P ,7(7p - %)) + tr - 7«)) + o(e). (4.17) 

To study the sign of (|4.17p . we remind that given an operator T and an orthogonal 
projector P, we can consider the block decomposition of T defined by 



T = 



FTP 


PT(1 - P) \ 


- i " + 


T+ 


(1 - P)TP 


(l-P)T(l-P) J 







(4.18) 



Moreover, let R be another orthogonal projector and consider Q = R — P; then 
P + Q is a projector and 



p + o 
p + o 

Q 2 



= {P + Qf 

= P + PQ + QP + Q 2 

= (l-P)Q-QP 

= (1-P)Q(1-P)-QP + (1-P)QP 

= Q—-Q++. 



(4.19) 



As a consequence, if Q = 0(e), then Q++ = 0(e 2 ) and Q-- = 0(e 2 ) 
Since 
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for fi = p,n, then we have 

+ tr (h~jA-~ (<£ - 7 P )A^) + tr (h^A+-(^ - 7n)A+ 7 
+ tr (h-^A-^1 - 7n)A" J + o(e) 

:= T p + + T- + T+ + T- + o(e) . (4.20) 

First of all, we analyze the relation between A^~, A^~ E and A^ 7E for /i = p,n. 
Using (13. lip , we obtain 



A ± =A ± 
A ± - =A ± . 



0(e), 
0(e), 
f O(e). 



Then, if we take P = A+ 7 e, = A+^-A+ 7 « and we apply KWi - KWi . we obtain 



A M.7 = K.Y- 



' 0(e 2 ) 


O(e) ' 


0(e) 


0(e 2 ) 



for /i = p,n. 

Moreover, since 7 e G jy, we can write 



0(e 2 ) 


0(e) 


0(e) 


1 + 0(e 2 ) 



7m++ 








o 



and 



V>7 'M V,7 



0(e 4 ) 


0(e 3 ) • 


0(e 3 ) 


0(e 2 ) 



So we can conclude that T p — o(e) and T n = o(e). 
Next, we remark that 

= tr (Kn A U% ~ %Kn) + tr (<7 A .,7^ ~ %Kn) • 

To calculate tr (-ff^ 7 A+ 7 (7^ — 7^)A+-Y we consider the block decomposition of 
A+- - A+- e for P = A+ - e . As before, we have 



A+ =A+ 



0(e 2 ) 


O(e) ~ 




0(e 2 ) 



l + 0(e 2 ) 


O(e) ~ 


0(e) 


0(e 2 ) 



for /i — p,n. 

Now, we observe that, in general, 7^ — 7^ = 0(e) and, more precisely, A^ ~ e (7^ 
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7/x)A* a = 0(e 2 ). Indeed, using the definitions from (|4.10|l to (|4.16|) . we have 



i 

Then 

M>7^ ifi I M> Mj7 

and 



0(5 2 ) 


0(£ 3 ) 


0(s 3 ) 


0( £ 4 ) 



^ + = tr (h+~A+~(% - 7m)A+ 7 ) + 
Next , we consider 7^ — 7^ . By definition, 

= (l - ^ [< 7 , 7/x] ) 7m (l + e + 7 , 7m] ) - % + o(e) 



Then 



= — e 



T+ 

M 



# m , 7 ,7m 
-etr (ff 



,7m 



(e). 



M>7 



^7'7M 



,7m 



,7m 



*t2 



+ o(e) 
+ o(e) 



(4.21) 



for [i = p,n and 

£(7p,7n)-£a,7«) - - £ E tr ( ff ^[[^'^ 

= 2e E ^ (« 7 7m) 2 " (K,) 2 

fi—p,n 

+o(e) 

where {A, B) = ii(A*B) is the Hilbert-Schmidt inner product. 
Then, using the Cauchy-Schwarz inequality, we obtain 

and 

£(7|>7n) -£(7p>7n) < °; 
furthermore, the equality holds if and only if (-ff^Vj^)* = iH^-j^. 
First, we consider the case (-H^Tm)* = H t-.^^ tllis im plies = tnat 



means 



#mVy'7m 



0. Then we have a contradiction. 
Second, if {H+-%)* = -#^ 7 7m, then 

7m^ 7 + ^7m = 
7m< 7 + %< 7 7m - 
7m^7 - = 
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7^ for /i = p,n. 

^ for ji — p,n, then we can construct 



that contradicts the hypothesis 

Finally, we can conclude that if 
7 e € Y\ n such that 

£(7p,7«) - £(l P ,ln) < 0, 
and thus we have a contradiction with the fact that 7 minimizes the energy on 



r 



Z.N- 



This implies that 



must be equal to zero and, as a consequence, 



for /i = p, n. 

As a conclusion, if g a -,gu-,9p and e are sufficiently small, "J is a solution of the 
equations (fl~26f and (|H?7|) . 

Appendix A. Proofs of lemma RQ1 and corollary 13.51 
In this section we give the proofs of lemma |3~H and corollary [ 



Proof of lemma \3^\ Given anMxM matrix B = (bij), we denote <& • B the right 
action of B on $ = (ipi, . . . ,<pm) G (L 2 (R 3 )) . More precisely, 



(M M \ 

i=l i=l / 



and, by straightforward calculation, we obtain 

Gram L2 ($ • B) = B*Gtaxa. L *{$)B 

where B* denotes the conjugate transpose of B. 
First of all, for [i = p,n, we consider 

and we observe that 



(A.l) 



Gram L 2 (# p ) 
Gram L 2 f# n J 



In- 



Second, we define 

$+ = 
p 



Gram i2 (A+ * p 



-1/2 



= A£**„« Gram L 2 A+ Vt *„ 



(A.2) 
(A.3) 



Remark that Grainy [Kp^pj and Grainy yA^^nj are invertible matrices 
thanks to the hypothesis [m] of the lemma. 

Next, we look for G (AT/^ii" 1 / 2 ^ x (a~^HV 2 ^ such that, taking 
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wc have 



0, 

0, 



(A.4) 
(A.5) 



with $ = ($ p , $„) and 1^,+ , defined by 



for /i = p,n. 



hi + Gram L2 + 



-1/2 



We observe that ^|,+ and are smooth maps from (A 9 H l ^\ to (i/ 1 / 2 ) and 



from ^A ?i ^i? 1 / 2 j to (i? 1 / 2 ) respectively; furthermore, 

Gram L 2 ($ p ) = G p , 
Gram L 2 ($„) = G„. 

Now, to prove the existence of <I>~ and $~, we apply the implicit function theorem. 
We remark that the equations (j A.4|) and (|A.5|) can be written as F(g, $~, <f>~) = 
where F is a nonlinear C 1 operator and g = {g^^g^^gp^e). In particular, 



and 



1 + Gram L 2 ( $ 



1 + Gram L 2 (jfr^ 

1/2 

• G} 2 . 



-1/2 



Hence, we define 



where 



-1/2 



F p ( 5 ,#-,$-) = $-. l + Gram L 2 ($-) • G p / 2 + X p (g, (A 



-1/2 



G^ 2 + (A.7) 



F„(<7, $ p , $„ ) = $„ • |1 + Gram L 2 ($„ 
and 

__ _ 1 f+°° 

K IM {g,<S> v ,$ n )=k^— J (H^-i^iH^-H^iH^-iTiy^^dr) 

for /i = p,n. 

Using the definitions (|1 .26[) and 271) . we obtain 

Xp(0, $ p , *-) = K n (0, = 0, 

and then F(0,0,0) = 0. 

Now, to apply the implicit function theorem, we have to check that 

z / , \ N 



" J • f-V,,// : (A, ,,// ; -' N V 



A ;,^ 1/2 ) x (a-^v2) j 
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is a C 1 operator and I?2-F(0, 0, 0) := F^- (0, 0, 0) is an isomorphism. We remark 
that 

( # q1/2 \ 
X " 1/2 ' ( A ' 8 ) 

1 /2 1 /2 

and then it is an isomorphism, since Gp'' and G n are invertible matrices. Pro- 
ceeding as above, we can easily show that F is well defined in ^A~ 1/<2 ^ x 

Next, we have to prove that F(g, <&~, $~) is C 1 ; by classical arguments, it is enough 

/ , \ z / , \ N 

to show that for ( X ,r) £ (A^tf 1 ^ x (A^if 1 ^ 

8% VP.* J - 

^ G i A -^ 



W p xe\S* H ') ' 

gMjjj rZ U- H vA N 



and we leave the details of this part to the reader. 

Then, applying the implicit function theorem, we conclude that there exist U C R 4 , 
V p C I A nil 1 ' 2 ) and V n C (A ' 2 ) neighborhoods of 0, and a unique 



continuously diffcrentiable function f : U ^ V p xV n such that F(g, f(g)) = 0; that 
means that for gu.g^.gp.e sufficiently small, there exists (i>~,$~) € V p x V n such 
that 

A p^ A p^p = °' 
K^K^n = 0. 

In particular, U = B(0, 7), VJ, = B(0,r]) and = B(0,r]) with 7,77 > and from 
the proof of the implicit function theorem, we know that, fixed 77, we can choose 
7 such that / : U — > Vp x V^. Then we take 77 and 7 such that D 2 F(g,x,T) is 
invertible V(#, x,t) G ?7 x V p x V„. 

Now, we denote J5 p := Gram L 2 ( A+^pj 



and we remark that 



*p = A ^**p + A p ,**p 

So we may write 

^.S-^^+A-^.S" 1 . 

As a consequence, 

^+( A p~,**p « V) = (*p ' ^ ' [ Gram i 2 (*p ' Bp')] 
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We can easily compute 

Gram L2 (* p . B~ x ) = (B^Graxa^^B- 1 = (BpB*)^ 1 
where B* denotes the conjugate transpose of B p . Since B p is hermitian, 
Gram L2 (^ p . S" 1 ) = (B^ 1 = (B; 1 ) 2 , 

k+(K,v*P • B p 1 ) = (*P • ^ « ( B p) 1/2 = *p 

«*+(A;,**p«Bp 1 ) -G! i /2 = *p- 



and 
Hence 



l*p-*pll(ffl/2)2 



with \I> p = A s \? p . In the same way, 



||*n-*».||(HV>)J» = 

with = A~<p*„ and B n := 



1/2 



( ff l/2)N 



(A.9) 



(A.10) 



Gram L 2 (A+ W * n ) 

We remind that the maps and are smooth; then, to have an estimation of 
the norms (|A.9|I and (|A.10[) , it is enough to estimate 



V V V (H1/2 y 

Indeed, Ve > 0, 3S p ,5 n > such that 



and 



n n n 



(JJV2)JV 



and 



P P P 



^ n ^ n 



(J?V2)2 



< <5„ 



(J? 1 /2)N 



<<*n 



(J?V2)2 



(H'/>)' 



< £ 



< e. 



Now, for 5 small enough, (^>~ «.B~ 1 , #~ •B" 1 ) eV p x V n ; then F(g, #~ »S~ 1 , *~ • 
B" 1 ) is differentiable and D 2 F(g, • B p • B" 1 ) := Q is invertible Vg £ £/. 

Using this fact, we can write 

F(g, = F(ff, • B- 1 , • B- 1 ) + Q($- - • B" 1 ) + u(g, S", 3") 

with 

*" = ($",*-) = /(<?), *" = B 

and 



B p 
S„ 



(if!/2)A 



lim 

and this implies 

Moreover, thanks to (|A.11|) . we know that there exists 5 > 0, such that 

1 



l«(s,2/)ll(ffi/3)/ 



< 



2IIQ- 1 ! 



y - * • B~ 
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if 



y- f " • B- 1 



(Hl/2-jA 



< 5. 



Then, choosing r\ < |, we have 



(if!/2)A 



< HQ- 1 !! F^^.fl- 1 ,*-.^- 1 ) 



( ff l/2)A 



(i/!/2)A 



and 



l>~ - • B^ 1 



< C 



< C5. 



(£fl/2)A 



Finally, choosing 5 < """(^j" 15 ") | we obtain 



and 



|$n - *n||(jyl/2 ) iv < £. 



To conclude the proof of the lemma, we have to show that 

A p,* : ImA p,$ "> ImA p,* 

A n,* : ImA n,* -> ImA n,* 

are one-to-one operators. We remark that 



A p,$ A p,* ^imA: 



A p,$ A p,* A p,$ 



A- A 



(A.12) 

(A.13) 
(A.14) 



p,* I il p>* p>* 



< 



A 



A p,* A p,$ 



< 



A p,* A p,* 



< 1. 



As a consequence, A p $ A p ^, is an invertible operator and A p ^ is one-to-one from 
ImA~ $ into ImA"^,. In the same way, we can prove that A~^, is one-to-one from 



ImA n,* int0 ImA n,*' 

In conclusion, 

A p,$$p - 0, 

K,**« = 0. 

□ 

Proof of coroUarv \3.5[ To prove this corollary, we apply lemma [3T4l to ^ fc for any 
k E N and, to obtain (|3.3ip and (|3.32p . we use the inequalities (|A.9|) . ()A.10|) and 

□ 
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